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ABSTRACT 

We derive the equations for the propagation of relativistic ionization fronts (I-fronts) in both static and moving 
gases, including the cosmologically-expanding intergalactic medium (IGM). We focus on the supersonic R-type 
phase that occurs right after a source turns on, and we compare the nonrelativistic and relativistic solutions for 
several important cases. Relativistic corrections can be significant up until the light-crossing time of the equilib- 
rium Stromgren sphere. When q, the ratio of this light-crossing time and the recombination time, exceeds unity, 
the time for the expanding I-front to reach the Stromgren radius is delayed by a factor of q. For a static medium, 
we obtain exact analytical solutions and apply them to the illustrative problems of an O star in a molecular cloud 
and a starburst in a high-redshift cosmological halo. Relativistic corrections can be important at early times when 
the H II regions are small, as well as at later times, if a density gradient causes the I-front to accelerate. For the 
cosmologically-expanding IGM, we derive an analytical solution in the case of a steady source and a constant 
clumping factor. Here relativistic corrections are significant for short-lived, highly-luminous sources like QSOs 
at the end of reionization (z ~ 6), but negligible for weaker or higher-redshift sources. Finally, we numerically 
calculate the evolution of relativistic I-fronts in the presence of small-scale structure and infall, for a large galaxy 
undergoing a starburst and a luminous, high-redshift QSO. For such strong and short-lived (z ~ 7) sources, the 
relativistic corrections are quite significant, and small-scale structure can decrease the size of the H II region by 
up to an additional ~ 25%. However, most of the IGM was ionized by smaller, higher redshift sources. Thus the 
impact of relativistic corrections on global reionization is small and can usually be neglected. 

Subject headings: H II regions — ISM: bubbles — ISM: galaxies: halos — galaxies: high-redshift — intergalactic 
medium — cosmology: theory 



1. INTRODUCTION 

Ionization fronts are ubiquitous in astrophysics. On interstel- 
lar scales, they occur when new stars form in molecular clouds, 
and their propagation has important dynamical effects on the 
surrounding medium. The compact and ultra-compact H II re- 
gions that form in this way have been extensively studied both 
theoretically and observationally (for recent detailed reviews 
and further references see e.g. Garay & Lizano 1999; Evans 
1999; Churchwell 2002). Radiative feedback effects play a key 
role in determining the physical conditions in star-forming re- 
gions, with profound effects on the star-formation rates and ef- 
ficiencies, and the resulting initial mass function. For example, 
ionization fronts may trigger star formation in dense clumps 
through compression and subsequent cooling of the gas (e.g. 
Lefloch et al. 1997; Hosokawa & Inutsuka 2005; Lee et al. 
2005). 

On larger scales, the propagation of ionization fronts has 
emerged as one of the key problems in cosmology. Shapiro 
& Giroux (1987) pointed out that the first sources of ioniz- 
ing radiation in the post-recombination expanding universe, be 
they quasars, primeval galaxies, or some pregalactic objects like 
Population III stars, would have caused cosmological I-fronts 
to sweep outward into the surrounding intergalactic medium 
(IGM), and they derived the equations to describe their propa- 
gation. These I-fronts, they found, were generally weak, R-type 
fronts that moved supersonically with respect to both the neutral 
gas ahead of them and the ionized gas behind, thereby racing 



ahead of the hydrodynamical response of the ionized gas. The 
H II regions created by these expanding cosmological I-fronts 
generally do not fill their "Stromgren spheres" - i.e. a fully- 
ionized sphere whose size is just large enough to encompass 
as many atoms recombining per second as there are ionizing 
photons emitted per second by the central source 4 . 

The origin of this difference between the cosmological I- 
fronts and their interstellar counterparts is the ever-decreasing 
mean gas density in the cosmological case. The predominance 
of this R-type phase in the cosmological case has made it pos- 
sible to approximate the propagation of the global I-fronts be- 
lieved responsible for the reionization of the universe by a 
"static limit" in which the evolving inhomogeneous density 
field of the IGM is assumed to be predetermined by large-scale 
structure formation, neglecting the hydrodynamical response 
of the gas to its ionization (e.g. Ciardi et al. 2000; Sokasian 
et al. 2003). On small scales, where gas can be much denser 
than average, such as inside halos, this "static" approxima- 
tion can break down as the I-fronts decelerate to about twice 
the sound speed of the ionized gas and transform from super- 
sonic R-type to subsonic D-type, thereby coupling their fate 
to the dynamical response of the ionized gas. Some attempt 
has also been made to study this dynamical phase of cosmic 
reionization I-fronts, by "zooming in" on the gas dynamics 
of the I-front as it encounters individual halos (e.g. Shapiro, 
Iliev, & Raga 2004; Iliev, Shapiro, & Raga 2005b) or by re- 
laxing the "static approximation" referred to above while sim- 
ulating large-scale structure formation during global reioniza- 
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tion (Gnedin 2000; Ricotti et al. 2002). Here the primary is- 
sue has been to characterize the onset, duration, and topol- 
ogy of reionization (Gnedin 2000; Bruscoli, Ferrara, & Scan- 
napieco 2002; Haiman & Holder 2003; Wyithe & Loeb 2003; 
Ciardi, Stoehr, & White 2003; Furlanetto, Zaldarriaga, & Hern- 
quist 2004; Iliev, Scannapieco, & Shapiro 2005a), particularly 
in light of the discrepancy between the epoch of reionization 
z <; 15 favored by the first-year data from the Wilkinson Mi- 
crowave Anisotropy Probe (Kogut et al. 2003), and observations 
of a Gunn-Peterson trough in high-redshift quasars, which sug- 
gest that reionization did not end until z ~ 6 (White et al. 2003; 
Fan et al. 2004). Other important questions include the growth 
of individual ionized regions around high-redshift Lyman-alpha 
emitters (Haiman 2002; Santos 2004), which would otherwise 
be highly obscured by neutral gas (Malhotra & Rhoads 2004; 
Stern et al. 2005), and the relationship between the proximity 
effect around high-redshift quasars and the mean IGM neutral 
fraction at this redshift (Wyithe & Loeb 2004; Oh & Furlanetto 
2005; Yu & Lu 2005). 

Until recently, the I-fronts described above have been treated 
nonrelativistically, as if the speed of light were infinite. Unlike 
shock fronts, which can reach relativistic (Rel) speeds only in 
the most extreme environments where the energy density is very 
high, such that the kinetic energy per baryon is comparable to or 
greater than its rest mass energy, ionization fronts are intrinsi- 
cally radiative phenomena that are mediated by the propagation 
of photons. This means that even the weakest ionizing source 
may, in principle, drive a strongly relativistic I-front. Thus the 
finite speed of light is always an issue for I-front propagation, 
whose importance should be ascertained for each set of astro- 
physical circumstances. 

Nevertheless, relativistic I-fronts have received only limited 
treatment in the literature. Roughly speaking, the I-front will 
be relativistic whenever the photon to atom ratio at the front 
exceeds unity. This occurs whenever a source releases ioniz- 
ing photons faster than a surface moving outward at the speed 
of light can overtake as many neutral atoms. The equations of 
Shapiro & Giroux (1987) for the evolution of the spherical I- 
front around a point source in the cosmologically-expanding 
IGM show that, at early times close to source turn-on, the 
nonrelativistic (NR) treatment formally yields an I-front speed 
which exceeds the speed of light, for just this reason. This "su- 
perluminal" phase is generally so short-lived, lasting only a tiny 
fraction of the age of the universe at source turn-on, that it has 
a negligible effect on the overall evolution of the I-front for 
a steady source. Recently, however, interest in this brief, early 
phase of the I-front has been drawn by observations of the spec- 
tra of high redshift quasars and Lyman alpha line emitters at 
redshifts z > 6 which suggest that we are seeing each of these 
sources as if it is surrounded by an isolated intergalactic H II 
region of its own making. In the quasar case, in particular, as- 
sumptions that the source of ionizing radiation is the observed 
quasar, itself, and that quasar lifetimes are much less than the 
age of the universe at those epochs, combine to suggest that the 
relativistic phase of the I-front predicted by the nonrelativistic 
solutions is quite relevant. For this reason, White et al. (2003) 
presented a revised equation for the expansion speed of the 
spherical I-front which takes account of the finite speed of light, 
although it neglects recombinations. They noted, however, that 
the relationship between the apparent size of the H II region 
along the line of sight and the parameters on which the physical 
evolution of that size depends, namely the quasar's luminosity 



and the average density of neutral H atoms in the IGM outside 
the H II region, are exactly the same as would be achieved by 
applying the equations of Shapiro & Giroux (1987). This sur- 
prising result is apparent also in the treatments of Yu (2005) 
and Yu & Lu (2005), who discussed the effect of relativistic 
expansion on the apparent shape of a cosmological H II region 
around a high-redshift QSO, including the effects of recombi- 
nation. It results from the fact that the photons in the observed 
quasar spectrum at the redshifted wavelength which marks the 
outer edge of the H II region were emitted by the QSO at the 
retarded time equal to the time those photons reached the I- 
front, minus their light-travel time from the source. Using this 
retarded time to evaluate the I-front radius in the equations of 
Shapiro & Giroux (1987) automatically accounts properly for 
the finite speed of light. 

These recent discussions of the QSO H II regions state their 
basic H II radius evolution equation without derivation and 
do not explain its connection to the fundamental conservation 
equations for an I-front, namely, the I-front continuity jump 
condition and the equation of radiative transfer for the ioniz- 
ing photons. In what follows, we will start from those fun- 
damental conservation equations to take full account of spe- 
cial relativity in generalizing them to address a wider range of 
circumstances and make clear the derivation that underlies the 
evolution equation for the spherically symmetric H II regions 
mentioned above. We will examine the general problem of rel- 
ativistic I-fronts, in both static, moving, and cosmologically- 
expanding media, deriving a number of exact and universal so- 
lutions and assessing the importance of relativistic effects in 
both the interstellar medium and cosmology. 

The structure of this work is as follows. The relativistic equa- 
tions for I-front evolution in a general form are derived in § 2 
for a source in a moving gas. In § 3 we use these equations 
to describe the relativistic evolution of a spherical I-front sur- 
rounding a point source in a static medium, deriving some exact 
analytic solutions to this problem, for a uniform gas density and 
a density that varies as a power-law in radius. For illustration, 
we apply these to a front expanding around an O-star or OB 
association in a molecular cloud and to a front inside a cos- 
mological halo. We also present a numerical solution for the 
I-front around a point source in a plane-stratified medium. In 
§ 4 we generalize the equations of Shapiro & Giroux (1987) for 
the cosmological I-front around a point source in an expand- 
ing IGM to take account of special relativity, finding an exact 
analytical solution, as well, which we use to assess the impor- 
tance of relativistic effects as a function of redshift and source 
flux. In § 5, we present numerical solutions for the isolated 
H II regions around cosmological sources, which include the 
effect of small-scale cosmological structure (i.e. self-shielding 
minihalos, time-varying clumping factor for the diffuse IGM, 
and the infall of gas surrounding each source halo) on these in- 
tergalactic I-fronts. We apply these solutions to determine the 
impact of relativistic effects on global cosmic reionization, by 
the approach described in Iliev et al. (2005a), which evolves the 
I-fronts created by a statistical distribution of cosmological ion- 
ization sources to the point at which their H II regions finally 
overlap. Our conclusions are summarized in § 6. 

2. BASIC EQUATIONS 

2.1. Relativistic I-front in a Moving Gas 

Consider the motion of an I-front through some gas that is it- 
self moving with respect to the source of the ionizing radiation. 
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In the lab frame defined as the rest frame of the source, the I- 
front moves with velocity V/ = /3/C, while the gas just outside 
the I-front on the neutral side (side "1") moves with velocity 
v g j = [3 g .\c. Assume that V/ and v gi i are in the same direc- 
tion, so the motion is one-dimensional. The relativistic law for 
the addition of velocities then gives the velocity of the I-front 
as measured in the rest frame of the gas just outside the front, 
v'j = (3jC, according to 



vi-v g .i 



l-v/v g ,i/c 



.2 " 



(1) 



In the rest frame of the I-front, the gas moves with velocity 



-f3'jC, where here and below double primed 
variables are quantities measured in the frame of the front. 

2.2. I-front Continuity Jump Condition 

The relativistic I-front continuity jump condition, which ex- 
presses the balance between the flux of ionizing photons arriv- 
ing at the front and the flux of atoms that are ionized by those 
photons as the atoms pass through the front, is written in the 
rest frame of the I-front as follows: 



(2) 



where v' g ' j is the front velocity relative to the undisturbed gas 
into which the front is moving, given by v'j [ in equation (1). 
To get the density (ijj,,we take the value measured in the rest 
frame of the I-front, 



(3) 



where the Lorentz factors, 7, are related to the velocities, (3c, 
by 7 = (1 -/3 2 )" 1 , as usual. The quantity /?,- counts the num- 
ber of ionizing photons required to yield a single newly ion- 
ized H atom on the ionized side of the front and is thus frame- 
independent. 

The ionizing photon flux F" as measured by the moving front 
is different from the lab-frame flux F measured in the frame of 
the source, by a Lorentz-boost to the frame moving with veloc- 
ity v/ in the lab frame: 

F"= 7/ (l-/3/)F (4) 

(Shapiro, Milgrom, & Rees 1986), where we have implicitly 
assumed that the relativistic Doppler shift of the photons from 
the lab frame to the rest frame of the moving gas does not sig- 
nificantly affect the integrated number of ionizing photons in 
the incident flux. Otherwise, we should replace F in equation 
(4) by the modified flux which counts only those photons at 
the location of the I-front at lab frame frequencies v such that 
v > Vmin, where 



V m in = ^(/,7g,l(l+/3 ? ,l), 



(5) 



where v is the threshold frequency for ionizing H atoms from 
the ground state. 

Combining equations (1) and (2)-(4) yields the equation for 
the I-front velocity in the lab frame as follows: 



(l-/?/&,i)7 g 



7/d -PiW 



(6) 



We have also used the fact that the lab frame H atom density 
n H j is related to the value in the rest frame of the gas, accord- 
ing to 

«ff,i=»H,i7s,i- (7) 



2.3. Radiative Transfer of Ionizing Flux Across H II Region 

2.3.1. Optical Depth 

Ionizing photons traveling radially outward from the central 
source will suffer absorptions due to the residual bound-free 
opacity of gas on the ionized side of the I-front, inside the H 
II region. As these photons cross a spherical shell of radius 
dr as measured in the lab frame (i.e. rest frame of source), 
in which frame the H atom number density is «#(?") and the 
frame-independent neutral fraction is (l—x), the probability that 
a photon emitted at a frequency v in the lab frame is absorbed is 
given for a shell that expands away from the source at velocity 
v g( r ) = 0g(f)c> by optical depth dr, according to: 

dr v = n H (l -x)a v t{\-(3 g )dr, (8) 
where ov is the photoionization cross section of an H atom in 
the ground state, evaluated at photon frequency v' = ^ g {\—(3 g )v, 
the frequency in the comoving frame of the expanding shell 
(Shapiro et al. 1986). 

2.3.2. Ionization Equilibrium 

For gas in the expanding shell, ionization equilibrium in 
the "on-the-spot" approximation, in which radiative recombi- 
nations to atomic levels n > 2 are balanced by photoionizations 
from the ground state, gives the neutral fraction according to 

l-x a B (n' H ) 2 a B n' H 



(9) 



x z n' H oF' oF< ' 
where a is the frequency-averaged bound-free cross-section for 
frequencies v' > v,h, weighted by the flux of photons at each 
frequency, and where the primed quantities are measured in the 
comoving frame of the gas. The atomic number density and the 
photon flux in this frame are related to their counterparts in the 
lab frame according to 

n> H = n H/lg, (10) 

and 

F' = lg (l-p g )F, (11) 

where F should be the photon flux for v > J g (l+(3 g )v t h (Shapiro 
et al. 1986). Equations (9)-(l 1) yield 
1 -x = a B n H 
x 2 <7 7 g 2 (l-/?g)F- 
According to equations (8) and (12), the optical depth through 
the shell is given, therefore, by 



2 2 

a B x n L H 



dr. 



(13) 



2.3.3. Flux Variation With Distance From Source 

Consider the lab frame quantity 5(r,f s ) = 4Trr 2 F(r,t), which 
is the rate at which ionizing photons pass through a shell of ra- 
dius r at time t that were emitted by the source at the retarded 
time tR = t — r/c. The quantity 5(0, f) is the ionizing photon 
luminosity of the source at time t. If the source luminosity is 
L v (ergs -1 Hz" 1 ) in the lab frame, then 

1 U 



S(0) = JV 7 



N~. v dv 



hv 



-dv. 



(14) 



In the absence of opacity due to photoionization, the quantity 
S(r, f«) is independent of r. The effect of opacity is to attenuate 
5(r,ffi) as follows, using equation (13): 



±[A^F v {r)]=-N^ v {t R )e-^ d -^ r 



= -[4nr 2 F v (r)] 



a B x 2 n 2 H 



■(15) 
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Multiplying both sides of equation (15) by dv and integrating 
over all frequencies v > 7 g (l +(i g )v t h yields 



^ = -4-Kr 2 a B x 2 n 2 H ^ g 



(16) 



We can integrate this over radius to yield S(rj,t R j) at the loca- 
tion of the I-front according to 



7" 

Jo 



S(n,t R j) = S(0,t RJ )- / dr-4Trr z xn e n H j g z a B , 



(17) 



where t R j = f«(r = r/) and we evaluate the quantities inside the 
integral on the right-hand-side at the time t(r) = t R j + r/c. We 
have also replaced one of the factors of xnn by n e above. Al- 
though we have neglected the presence of helium until now, 
we can account for the contributions of ionized helium to the 
electron density inside the H II region by writing n e = Xe«xnn, 
where XeS = 1 + pA(He), A(He) = 0.08 is the He abundance by 
number with respect to hydrogen, and p = 0, 1 or 2 according to 
whether He is mostly neutral, singly ionized or doubly ionized, 
respectively. If we also allow for a volume-averaged clumping 
factor, C = (riff) / (hh) 2 , then equation (17) becomes 



S{n,t R ,)=NJt R ,)- 



dr-4irr 2 XeffX 2 n H j 2 a B C. (18) 



Here a B = 2.59 x 10" 13 cm 3 s" 1 is the case B recombination co- 
efficient for hydrogen at 10 4 K (which is the temperature we 
assume for the ionized medium throughout this paper). The 
value of the radius rj which makes S(rj,t R j) vanish according 
to equation (18) is then the Stromgren radius r s . Hereafter, un- 
less otherwise noted, the only retarded time to which we shall 
refer will be t RJ , and we shall drop the subscript "I" when refer- 
ring to it. The lab frame flux at the front at time t is then given 
by 

S(r,,t R ) 

F(n,t)^ 



4nrf 



(19) 



While the analysis in § 2 apparently applies to a point source 
of ionizing radiation surrounded by a spherically-symmetric 
gas, it actually applies equally well to the general case of in- 
homogeneous gas which varies in all three spatial dimensions. 
In the latter case, we simply interpret the equations as if they are 
written in (r, 9, (f>) spherical coordinates centered on the source, 
as long as any non-radial gas motions which might be present 
are subrelativistic. If a problem arises in which relativistic, non- 
radial gas motions occur, the equations above should be appro- 
priately modified to account for the effect of transverse velocity 
on the relativistic transformations, as well. In the absence of 
such relativistic transverse gas motion, however, we can use the 
analysis presented in this section for the general 3D problem 
without modification. 

3. RELATIVISTIC I-FRONTS IN STATIC MEDIA 

3.1. Point Source in a Static Medium 

Consider a point source in a static gas. In this case, v gi i = 0, 
the velocity of the I-front v\ in the frame of the undisturbed 
gas ahead of the front is just the same as its velocity v/ in the 
lab frame in which the source is at rest, and 7^ = 7/. In the 
frame of the front, therefore, v" g X = -V/, the H atom density 
n H 1 = li n H,\, and the flux F" is related to the lab frame flux 
F by equation (4). Equation (6) for the front velocity in the lab 
frame then becomes 

fiin HA c+F 



where F is the number flux of ionizing photons at the current 
position of the I-front, given by equation (18) with 7 g = 1 and 
equation (19). 

Combining equations (18) - (20) and solving for v 7 , we ob- 
tain the lab frame evolution equation for the radius r; of the 
spherical I-front, 



dn _ N 7 (t R )-4ir Jq' n H 2 Ca B x e &r 2 dr 



(21) 



4nf3iCrjn H ^ +Ny(t R )-4ir J Q ' n H 2 Ca B x e «r 2 dr 

In the nonrelativistic limit in which dr; jdt <C c, the first term in 
the denominator becomes dominant and equation (21) reduces 
to the usual expression 



dr; 
dt 



A> 7 (0-4tt J f ' n H 2 Ca B x e sr 2 dr 
4n(3irjn HA 



(22) 



In the opposite, strongly relativistic limit, in which the flux is 
very large, both the denominator and numerator are dominated 
by the N y term and dn jdt approaches c asymptotically. 

For this static medium, the I-front velocity equation can also 
be derived by imposing photon conservation within the ionized 
volume. This is done by requiring that at a time t the total num- 
ber of photons that have reached the ionization front is equal to 
the sum of the total number of atoms within the volume and the 
total number of recombinations that occurred within the volume 
excluding recombinations that affect emitted photons that have 
not yet reached r/(f ), that is 



/ N 7 (t')dt' = 47T/3,- / n H (r)r 2 dr 
Jo Jo 

fdt'f 

Jo Jo 

I 



+ 4-7T 



n 2 H (f)Ca B x&r 2 dr 



■4tt/ dt' I n 2 H (r)Ca B x eS r 2 dr (23) 
10 Jo 

Differentiating this equation with respect to time, we obtain 



N 1 (t R )-4ir I n H (r)Ca B x e ffr dr 




n 



drj 



(24) 



and solving for fj recovers equation (21). 

For a given density field and time-dependent photon lumi- 
nosity, equation (21) can, in principle, be solved to yield the 
I-front radius and radial velocity as explicit functions of time 
t, for any direction (0,<fr). We will do so for certain density 
profiles in the sections that follow. Before we consider such ex- 
plicit solutions, however, we note that there is a simplification 
that results if we rewrite our equations in terms of the retarded 
time, t R , rather than physical time, t. For every interval dt of 
physical time, there is a corresponding interval 

dt R =dt(l-p I ). (25) 

If we define a "retarded velocity" /3i,r = c~ l drj/dt R , then equa- 
tion (25) says that /?i jR is related to (3j according to 



(26) 



l-AW 

Combining equations (19), (20) and (26) then yields an equa- 
tion for dri/dt R in terms only of t R , as follows: 

dn S(n,t R ) 



dt R 4TTrjn H ,iPi' 



(27) 
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Remarkably enough, equation (27) is identical to the nonrela- 
tivistic equation for the I-front velocity drj/dt, equation (22), 
if we replace t in the latter equation by t R . 5 Therefore, any 
solution of equation (27) for the explicit dependence of r/, and, 
using equation (26), of /?/, on t R becomes an implicit solution 
for the dependence of rj and (3i on physical time t by relating t 
to t R according to 

t = t R + r,(t R )/c. (28) 

This universal property will allow us to simplify the solution of 
the I-front evolution equation and yield analytical results when 
the explicit equation (21) otherwise would seem to require a 
numerical solution. 

If we specialize to the case in which the source luminosity 
after turn-on, 7Y 7 , is time-independent, then, according to equa- 
tions (18) and (27), the "retarded velocity" f3i. R , when the I- 
front radius has some value r/ is identical to the velocity vi.nr 
of the nonrelativistic I-front solution at the same radius, divided 
by c. In that case, we can write the dependence of the relativis- 
tic I-front velocity v/ = cj3; on r/, using equations (26) and (27), 
directly in terms only of the nonrelativistic solution at the same 
radius, according to 

= t~i — — -rvr ' < 29 ) 

l+V/^(r/)/c 

regardless of the difference in physical times at which the I- 
fronts in the relativistic and the nonrelativistic solutions achieve 
the radius r t . 

The relativistic I-front solution for the radius versus time, 
which depends on the source luminosity and the atomic density 
in the external medium, describes the size of the H II region 
in the source frame. In order to use this solution to interpret an 
observed H II region, however, we must take account of the fact 
that photons which travel from the I-front to the observer have 
different travel times if they originate from different locations 
along the front surface but arrive at the same time. Photons 
received from the part of the I-front which moves towards us 
along the line of sight to the central source travel a smaller dis- 
tance to reach us than did photons arriving at the same time 
from the part of the I-front which moves transverse to that line 
of sight. As a result, the line-of-sight radius will appear to be 
bigger than the transverse radius, since we see it as it was, later 
in the time evolution of the spherical front. Suppose we use the 
apparent transverse size measured at the distance of the source 
to fix the time t (measured from source turn-on) at which the 
relativistic I-front solution matches this radius, r/(f). In that 
case, the apparent line-of-sight radius will be given, instead, 
by r IyNR (t), the solution of the nonrelativistic I-front equations, 
evaluated at the same time t, since the time t is then equal to 
the retarded time for the apparent radius of the I-front along the 
line of sight. This apparent departure of the shape of the I-front 
from spherical was also discussed in the context of quasar H II 
regions at high redshift by Yu (2005). A related argument also 
explains why the nonrelativistic solution for the I-front radius is 
valid for interpreting the apparent line-of-sight sizes of the H II 
regions around high redshift quasars based upon the absorption 
spectra of these quasars, as mentioned in our introduction. 

3.2. Steady Source in a Uniform Density Gas 

If we assume a steady ionizing source that turns on in a gas 
with uniform density, clumping factor, and Xeff, all indepen- 
dent of time, then equation (21) has a closed-form analytical 

5 A similar correspondence was noted by Yu (2005) (see equation [7] there). 



solution, which is obtained as follows. Let us define the dimen- 
sionless ionized volume as 

y = V,/V s = V,Ca B n 2 H x,ffN- x , (30) 

where V s is the the Stromgren volume, and the dimensionless 
time as 

w = t/trec = t Ca B n H x e s, (31) 

where t rec is the recombination time. The non-dimensional form 
of (21) is thus given by 

dy (1-y) 



dw l+q(l-y)/(3y 2 ^y 



where 



Ctrec 



(32) 



(33) 



and (4-n/3)rj = Vs, i.e. q is the light-crossing time of the Strom- 
gren radius in units of the recombination time. This equation 
reduces to the nonrelativistic one when c — > oo, q — > 0, 

dy 



dw 



= i-y, 



(34) 



which has the standard solution 

y=l-e~ w . (35) 

In general, the solution of equation (32) is given by the implicit 
relation 

w = qy l/3 -ln(l-y). (36) 

The reader can easily confirm that this relativistic I-front so- 
lution in equation (36) is identical to that which results if we 
replace t in equation (35), the solution of the nonrelativistic I- 
front equations, by the retarded time t R , as described in § 3.1. 

In Figure 1 we show the I-front radius r t and velocity v/ ver- 
sus time t for different values of q. When q = 0, we recover the 
NR solution, in which the I-front starts with infinite velocity 
and decelerates to zero velocity within of order the H II region 
recombination time, by which time the I-front radius has grown 
to equal the radius of the Stromgren sphere. When q > 0, the 
speed of the I-front remains close to the speed of light from the 
time of source turn-on to roughly a time qt rec later. For q > 1 , 
the I-front does not reach the size of the Stromgren sphere until 
a time qt rec after the source turn-on. In fiducial units, q is given 
by 



q = 0.019C 2 ^N l %n Hr 



1/3 



(37) 



(assuming Xeff = 1-08, i.e. He is once-ionized) where N^^g = 
N 7 /(l0 49 s~ l ) and n#,3 = n#/(10 3 cm~ 3 ). This estimate shows 
that the astrophysically-relevant values of q can be significantly 
different from 0. Hence, the relativistic correction to the stan- 
dard I-front propagation solution can be important, especially 
in the immediate vicinity of the source. The correction is even 
more important when the ionizing source is very luminous or 
the surrounding gas is dense and highly clumped. In such 
strongly relativistic cases, the correction would be important 
during the entire R-type phase of the evolution of the H II re- 
gion, until the corresponding Stromgren sphere is reached (see 
Figure 1). 

To estimate the impact of our results for compact H II re- 
gions in more detail, we consider an O-type star in a molecular 
cloud environment. Here the typical gas number densities are 
nil ~ 10 2 - 10 3 cm~ 3 , but can reach <~ 10 4 - 10 5 cm~ 3 in denser 
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FIG. 1. — Relativistic I-front for a steady source in a static, uniform gas: (a) (top) radius (in units of Stromgren radius r$) and (b) (bottom) velocity (in units of 
rsArec)- Curves are labeled by values of the dimensionless light-crossing time of the Stromgren radius, q = rs/(rt re c), with q = (i.e. nonrelativistic limit) and 
q = 0.01,0.1, 1,10, and 100. In these dimensionless units, the speed of light is q~^ . 
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regions and up to Hh ~ 10 7 cm 3 in so-called "hot cores" (Garay 
& Lizano 1999; Evans 1999; Churchwell 2002). Furthermore, 
the observed gas distribution is very inhomogeneous (Garay & 
Lizano 1999; Evans 1999; Lebron et al. 2001), and the fill- 
ing factors of dense clumps are / v = 0.03-0.3 (Evans 1999), 
roughly corresponding to clumping factors C « f~ l = 3-30. 
Thus, in our estimates below we assume the mid-range value 
of C = 10. Finally, for a typical massive O-type star the ioniz- 
ing photon flux is N y ~ 10 49 s _1 , and up to a few times 10 50 s _1 
for a cluster of such stars with corresponding H II region radii 
ranging roughly from 0.01 pc to ~ 1 pc (Garay & Lizano 1999). 

How important are relativistic effects in such an environ- 
ment? ForC= 10, - 10 4 V, and n H ~ 10 5 cm 3 , we obtain 
q = 0.41, while using iV 7 ~ 3 x 10 50 s _1 , and n# ~ 10 4 cm" 3 re- 
sults in q = 0.59. As seen in Figure 1, this means that the I-front 
expansion is relativistic all the way out to the Stromgren ra- 
dius, which means that the resulting H II region size during this 
rapidly-expanding R-type phase is much smaller than the one 
predicted by the nonrelativistic expression. 

On the other hand, this rapid expansion phase is com- 
pleted within a few recombination times or less, where f rec = 
1 /(CasnnXes) ^ 1 vr > which is quite short compared to the dy- 
namical time of such compact H II regions, fd yn ~ 5 x 10 3 yr 
(Garay & Lizano 1999). Hence most of their evolution is spent 
in the slow D-type phase, in which the hot, over-pressured ion- 
ized bubble expands into the surrounding medium, following 
the initial, supersonic R-type phase. During this longer-lived 
D-type phase, the I-front moves with speed V/ -C c, and the rel- 
ativistic corrections are negligible. Numerically, in terms of the 
H II region final equilibrium size, <J 5% of the expansion occurs 
during the initial relativistic phase. 

Our approach neglects the finite rise time of the source lumi- 
nosity. Yet, for the high densities typical of star-forming regions 
in molecular clouds, this is much longer than the recombination 
time. For a 30 M Q star, for example, model star calculations 
find that the UV luminosity beyond the Lyman limit doubles in 
4000 years as the star approaches the main sequence (Mathews 
& O'Dell 1969). In this case, there may not be a distinct rela- 
tivistic R-type phase as estimated above. If the source luminos- 
ity Nj increases gradually over a time which is long compared 
with the recombination time, then the dimensionless parameter 
q defined above must also increase gradually over time, in pro- 
portion to TV-J/ 3 , as does the instantaneous Stromgren radius. As 
such, the criterion that q is significantly larger than zero, neces- 
sary for relativistic corrections to be important during the entire 
R-type phase, may not occur before the I-front fills this gradu- 
ally increasing Stromgren radius. On the other hand, realistic 
density profiles for the immediate surroundings of a massive 
star when it turns on are likely to be far from uniform. If there 
is a steep enough negative density gradient outside the star, then 
the I-front can accelerate outward even if it initially deceler- 
ated and entered the D-type phase (Franco, Tenorio-Tagle, and 
Bodenheimer 1990). This phenomenon is associated with the 
"champagne phase" of the H II region. When this happens, the 
I-front moves ahead of the dynamical expansion of the inner 
H II region as a weak R-type front and can accelerate to rel- 
ativistic speeds again. In that case, the finite rise time of the 
stellar luminosity is not relevant, since the flux at the position 
of the I-front effectively rises in the time it takes light to cross 
a few absorption mean free paths in the neutral gas outside the 
front. We will discuss the propagation of relativistic I-fronts in 
density gradients in §3.3 and §3.4 below. 



3.3. Steady-Source in a Power-law Density Profile 

Consider a source that switches on in the center of a 
spherically-symmetric density profile n(r). The case in which 
the density decreases with increasing radius is relevant when- 
ever the source, like a massive star, is forming by gravitational 
instability, in the middle of a centrally-concentrated gas pro- 
file. As in § 3.2, we shall assume that the source luminosity 
is time-independent and that the clumping factor is constant in 
space and time. The nonrelativistic problem of an H II region 
in a spherically-symmetric density distribution that varies as a 
power of radius outside of a flat-density core was discussed by 
Franco et al. (1990). Let the H atom number density in the 
undisturbed gas be defined by 

"7 

(38) 



r > r , 
r < r . 



As long as the I-front is inside the core, its propagation fol- 
lows the solution for a uniform-density gas derived in § 3.2 
above. In this phase the I-front continually slows down from 
V/ « c at small radii rj. If the core Stromgren radius r s .o < ro, 
where rs.o = [3^V 7 /(47rCaB«g)] 1/3 < ro, then the I-front will 
slow down to zero velocity just as it fills this Stromgren sphere, 
thereby remaining trapped within the core. If r s .o > ro, instead, 
then the I-front continues to expand beyond the core radius. We 
shall assume from now on that rs.o > ro- 

The time evolution of the I-front radius is given by insert- 
ing the power-law density profile in equation (38) into equa- 
tion (21). We nondimensionalize the resulting velocity equa- 
tion by expressing all radii in units of rs.o and time in units of 
f rec ,o = (agnoC)' 1 , the recombination time in the core 

Y=—, (39) 
rs,o 

so Yj = rj/r S fl, Y = r /r S fl, and Y s = r s /r S fl, while 

t 



w = 



frec.O 



(40) 



We also define the dimensionless ratio of the light crossing and 
the recombination time in the core as 

(41) 

The solutions for the I-front radius and velocity for each value 
of the density profile slope 7 are then fully characterized by the 
two dimensionless parameters q and Yq, as follows: 



dY, 
dw 



27 ' 
.27-3, 



3 ■ 
.27-3, 



y2 7y 3-2 7 
1 t\ 1 1 



/7 v 2-7 



+q 



2 7 -3 

31n(7//y ) 



2 7 -3 



27 v 3- 2 7 



i)(^)%[^ 3 -l-31nW/yo)] ! 



for 7^3/2, 

for 7 = 3/2. 
(42) 

For comparison, the speed of light in these dimensionless units 
is q~ l . 

The I-front leaves the core with initial velocity 
dYj_ \-Yl 

y i=Yo 3Y Q 2 + q(l-Y^y 
which is relativistic (roughly) if q > (3F 2 )/(1 - F 3 ) = 

[vi, N R(?"o)/(?"s,o/frec,o)]" 1 , where vi iNR (r ) is the nonrelativistic 
I-front solution speed when rj = ro (i.e. take q=0 in eq. (43) 
above). Hence, for any given ro and no, a sufficiently high lu- 
minosity N 7 is required to make the I-front velocity still rel- 
ativistic once the front reaches n = ro- It is possible for the 



(43) 
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FIG. 2. — Relativistic I-front for a steady source in a static gas with a power-law density profile, n# oc r~T, and a constant density core for r < rg (left panels for 
7 = 2, right panels for 7 = 2.5): (a) (top) radius (in units of core Stromgren radius rs.o) and (b) (bottom) velocity (in units of rs,o/frec,o)- Curves are labeled by values 
of the dimensionless light-crossing time of the core Stromgren radius, expressed in units of the core recombination time, q = rs.o/(rfrec,o). In these units, the speed 
of light is q~ l . The case q = corresponds to the nonrelativistic solution. All curves assume Yg = 4 -1 / 3 . 
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I-front to accelerate afterwards, however, depending upon the 
values of 7 and To, so even if the front is not relativistic when it 
leaves the core, it may become relativistic at larger radii. 

Franco et al. (1990) derived some of the properties of the R- 
type I-front phase for the nonrelativistic solution of this prob- 
lem. As discussed in § 3.1, we can use this nonrelativistic so- 
lution directly to derive additional properties of the relativistic 
solution, as well. The nonrelativistic I-front has a velocity vlnr 
which depends upon its radius rj as follows: 



vi.nrOo) , . 
vi.nrO/) = —5 — —u(V, 



y-3_i 



where u(-f) is given by 



w(7) = 



y-3. 



(y 7 /y ) 1/2 [Y ( f-l-3lnY,/Y Q 



2^=3 + 2^3 ( Y l/ Y o) 3 ^ 



(44) 



.**7*3/ &) 
for 7 = 3/2. 



Equations (44) and (45) then yield the correct relativistic veloc- 
ity v/(r/) for the I-front if we insert this vi NR given above into 
equation (29). 

For any I-front that expands beyond the radius of the core, 
the relativistic I-front will only expand until it reaches the same 
Stromgren sphere radius r s as in the nonrelativistic solution, 
given by 



rs,o 



3-27 27, 

3 T 3 1 



2 ■ 



yoex P {i[y - 3 -i]}, 



for 7^3/2, 
for 7 = 3/2 



(46) 



(Franco et al. 1990). Finally, there is a (flux-dependent) critical 
value of the logarithmic slope of the density profile, -7/, below 
which there is no finite Stromgren radius rs at which the I-front 
is trapped, given by 



7/ : 



2 L 



(47) 



(Franco et al. 1990). For density profiles that decline more 
steeply than r~ 7/ , the relativistic I-front expands without bound, 
just as it does in the nonrelativistic solution. 

The relativistic and nonrelativistic I-front propagation solu- 
tions for power-law density profiles are plotted in Figure 2 for 
the illustrative cases of 7 = 2 and 2.5, for the particular value of 
Yq = 4~ 1//3 for which the critical slope 7/ = 2. In that case, for 
7 = 2, the I-front is never trapped at a finite Stromgren radius, 
but it decelerates continuously and reaches zero velocity at infi- 
nite radius. Such I-fronts are therefore relativistic only at early 
times. For 7 = 2.5 > 7/, on the other hand, the I-front reaches a 
minimum velocity and thereafter accelerates, so it is relativistic 
both at early and late times. 

For example, let us consider a source producing f 1 = 250 ion- 
izing photons per atom within a lifetime of t s = 3 Myr, inside a 
galactic halo of total mass 10 n M Q , switching on at z = 7. This 
corresponds to one of the illustrative cases for the cosmolog- 
ical H II regions of individual sources during the reionization 
epoch considered in Iliev et al. (2005a), except that here we are 
interested in the early phase of I-front propagation inside the 
halo, while Iliev et al. (2005a) were interested in its subsequent 
propagation in the external IGM. For cosmological halos in the 
current ACDM universe, the typical density profile logarithmic 
slopes outside the central region are steep, in the range -2.1 to 
-2.4 (e.g. Truncated Isothermal Spheres [TIS], Navarro, Frenk 
and White [NFW], or Hernquist profiles). Let us, therefore, as- 
sume 7 = 2.3 as a typical value, but, for illustrative purposes, 
adopt a core density no and a radius ro related to the halo virial 



radius r v j r as in the TIS model of Shapiro et al. (1999) and Iliev 
& Shapiro (2001), for which r v j r « 30ro. For this mass halo at 
this epoch in ACDM, ro = 0.7 kpc, no = 1.64cm" 3 , and we ob- 
tain A^ 7 = 4.616 x 10 55 s _1 and r s ,o = 8.1 kpc, implying q = 0.36 
and Yq = 0.086. The I-front velocities are r/ = ro and r/ = r v ; r are 
then v/(ro) = 0.941c and v/(r v i r ) = 0.978c, respectively. Hence, 
in this example, the I-front is relativistic inside the galaxy at 
all times and is accelerating once it exits the core. Even assum- 
ing an order of magnitude smaller rate of production of ionizing 
photons, the velocity is still relativistic, starting at v/(ro) = 0.62c 
as the front exits the core, and accelerating to v/(r v j r ) = 0.82c. 
Of course, the actual density profile of gas in such a halo may 
be affected by radiative cooling and by departures from spher- 
ical symmetry which are not accounted for in this estimate, so 
these numbers are only intended to be illustrative. 

3.4. Plane-stratified Medium 

Given that sources tend to form in regions in which gravita- 
tional instability leads to collapse, the medium may not only be 
centrally concentrated but may also depart from spherical sym- 
metry. Gravitational collapse tends to result in flattened struc- 
tures, either by a "pancake" instability or by the formation of 
a rotationally-supported disk. It is therefore instructive to con- 
sider the I-front for a point source in a static, plane-stratified 
medium, therefore. This problem was also considered in the 
nonrelativistic limit by Franco et al. (1990). 

The density profile is given by 



«(z) = « sech 2 (z/z ), 



(48) 



where no is the central density, z is the height above the central 
plane at z = 0, and zo is the scale height. Once again, we define 
the dimensionless parameters 

To = — (49) 

rs,o 

and 

rs,o 



CfrecO ' 



(50) 



where rs.o and f rec ,o are the central Stromgren radius and recom- 
bination time, respectively. Application of equation (21) to the 
stratified density profile of equation (48) leads to the following 
angle-dependent differential equation for the evolution of the 
I-front radius, 

dYj 1-/07,0) 



dw 3Y?g(Y I ,6) + q(l-f(Y I ,6)y 
where Y0) = r/(0)/r s ,o, w = t/t recfl , 



f(Y,9)- 



0.66T 3 
sin 3 (6») 



■ tanlr 



Tsin(6») 
0.88T 



(51) 



(52) 



and 



g(T,0) = sech 2 [7/y o sin(0)]. 



(53) 

In deriving equations (5 1)-(53), we have measured the polar an- 
gle 9 of the direction from the central source to a point on the 
I-front with respect to the z = plane, so 6 = tt/2 corresponds 
to the symmetry axis (i.e. the z-axis). We have also followed 
Franco et al. (1990) in evaluating the recombination integral in 
equation (21) along each direction by approximating the inte- 
gral using 

r " p 2 sech 4 (p)dp ~ 0.22tanh 3 ( JL) , (54) 



/' 

Jo 



where p = (r/zo)sin#. 
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FIG. 3. — Two-dimensional, axisymmetric I-front surfaces (xi /rs.o,zi /rs.o) (in units of the Stromgren radius at the central density) for a steady source in a static, 
plane-stratified medium at different times after turn-on, f/f rec .o (where ? rec o is recombination time at the origin), as labelled, for q = (i.e. nonrelativistic solution) 
(left panels) and q = 1 (relativistic) (right panels). Both curves assume Yq = 1/2. 
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FIG. 4. — Relativistic I-front radius (upper panel) (in units of Stromgren radius rsfl at central density) and velocity (lower panel) (in units of rs,o/frec.o) along the 
symmetry axis (z-axis) versus time (in units of central recombination time f rec ,o) f° r sam e plane-stratified case as shown in Figure 3. Each curve is for a different 
value of q = rs.o/(cf re c.o) (as labelled in lower panel), with q = corresponding to the nonrelativistic solution (solid). In these units, the speed of light is q~ [ . Vertical 
dotted line marks the finite time toe at which the NR solution (q=0) reaches infinite radius: zi(too) = oo. All curves assume Y = 1 /2. 
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Illustrative solutions of equations (51)-(53) for relativistic 
and nonrelativistic I-fronts in a plane-stratified medium are 
plotted in Figures 3 and 4. We adopt the value Yq = 1/2 in 
all cases. Figure 3 shows the two-dimensional, axisymmetric I- 
front surfaces at different times for the nonrelativistic solution 
(q = 0) and for the relativistic solution for q = 1 . The NR I- 
front solution starts out at superluminal speeds and decelerates 
in all directions. Along the symmetry axis, the NR I-front even- 
tually reaches a minimum velocity and, thereafter, accelerates 
upward. This acceleration leads to a superluminal "blow out" 
in which the NR I-front reaches infinite height in a finite time 
too. ln the perpendicular direction, however, along the plane 
of symmetry at z = 0, the NR I-front decelerates continuously 
and approaches the Stromgren radius of a uniform sphere of 
the same density. The relativistic I-front also starts out decel- 
erating but remains close to the speed of light in all directions, 
so its shape is initially quite spherical. Like the NR I-front, 
the relativistic I-front also reaches a minimum speed along the 
symmetry axis, before it accelerates upward once again and ap- 
proaches the speed of light. Since its speed is always finite, 
however, the relativistic I-front cannot "blow-out" as the NR 
front does in a finite time. Since, in the central plane, the rela- 
tivistic I-front slows to approach the same Stromgren radius as 
does the NR solution, above this plane it must balloon upward 
and outward, confined at the waist by this "Stromgren belt". 

4. COSMOLOGICAL RELATIVISTIC I-FRONTS 

Next we generalize the treatment by Shapiro & Giroux 
(1987) of an I-front in a uniformly expanding cosmological 
background at a redshift z, to take account of special relativ- 
ity. Here it is convenient to define the cosmological scale fac- 
tor as a = (1 +z,)/(l +z), where z,- is the redshift of the source 
turn-on, corresponding to initial time f,. The comoving hy- 
drogen number density («h,<)> I-front radius (r/), and veloc- 
ity (v/)> then correspond to quantities in proper coordinates 
(in the frame of the source) as n H)P = n H ^/a 3 , r^ p = arj, and 
V/j, = dr^p/dt = d(ari)/dt, which means that these comoving 
quantities are defined at the source turn-on time, not at the 
present as is often done. 

4. 1 . Relativistic I-front Velocity in the Cosmologically 
Expanding IGM 

For a relativistic I-front in the cosmologically expanding 
IGM at redshift z, the velocities of the front and of the gas just 
outside it as measured in the lab frame (i.e. the source rest 
frame) are the proper velocities, v Lp = fa tP c and v g j, p = (i g .\, p c, 
respectively. This latter gas proper velocity is the sum of the 
proper Hubble velocity at the location of the front at the proper 
radius r^ p , Vh(ri, p ) = fa(r^ p )c, and its proper peculiar velocity, 

Vg,pec> 

v g ,i,p = Vh(ri, p ) + v g)Pec 

= H(z)r l!P + v g . pecl (55) 

where H(z) is the Hubble constant at redshift z. We are inter- 
ested in I-front radii that are small compared with the horizon, 
so Hr I p <C c, and in situations in which v g , pec <cas well. For 
the problem of an I-front in the uniformly-expanding, mean 
IGM, in fact, we can set v g _ pec = 0. The I-front velocity mea- 
sured by the gas at that location is the proper peculiar velocity 
v' I pec , given by equation (1) as 



v I,pec ' 



vi,p-Vh(n, P ) _ Vi, p -H(z)r hp 

Vh(ri,p)vi.i, ~ ^ gfe)r/, p v/,p 



1 



(56) 



Let us define the proper peculiar velocity of the I-front as mea- 
sured in the lab frame of the source as V/ iPec = fa, pec c, where 

dr; 

vi. pe c = vi, p -H{z)r Lp = a—. (57) 



dt 



Equations (56) and (57) then yield 

^l.pec ~ fil.perC ~ 



I.pec^ 



i-A/V (58) 

According to equation (6), the proper peculiar velocity of the 
I-front in the lab frame of the source is given in terms of these 
quantities by 



I, pec : 



d-fa, P fa) 



'I,pec 

where fa tP = p\, ptc + fa, 7i' 



1iA 1 -Pi,p)F 



fan HtP c 



(59) 



_ 7(/3 I ' pec ), and we have replaced 
nu,\ by «H,p7/! (using equation [7]). Substituting equation (58) 
into equation (59) yields an equation that can be solved for p\ tPec 
as a function of fa and the quantity F /(/3in HtP c). There are two 
roots to this equation, corresponding to /3i )Pec > 1 and /3i pec < 1, 
respectively. The physical solution must have /3[ jPec < 1, which 
yields 

_adr, (l-0h)F 
Pi.pec ---77 = ~6 tt; • (60) 

c dt PiitH^hC+F 

The photon flux F in equation (60) is that at the location of the 
I-front as measured in the frame of the source, given by equa- 
tions (18) and (19) as 

S(n,t R ) 



F(r,,t): 



4-na 2 r 



(61) 



with 

S(ri,t RtI )=N 7 (t RjI )- 



dr p -4Trr 2 x e ffX 2 n 2 H ~/ h 2 a B C, (62) 



where the proper radius r p is related to the comoving radius r 
by r p = ar, the retarded time is f R i = t — i\ }V /c and the quantities 
inside the integral are evaluated at the time t' = t^\ + r p /c. 

The Hubble velocities across the H II regions of interest to 
us here are always much smaller than the speed of light. For a 
flat universe in the matter-dominated era, with total matter den- 
sity Clo m units of the critical density and Hubble constant h 
today in units of 100kms _1 Mpc _1 , cosmic expansion results in 
a Hubble velocity at the proper radius, r\ p , of the I-front, given 
by 

&( ^ ) = ^.4.MX,0-^VV+A , '7 



0.15 



1+z 
7 



\6Mpc 
(63) 

Since the largest H II regions (e.g. around rare luminous QSOs) 
are estimated to be several Mpc in proper radius at z ~ 6, this 
implies that f3/ t is always less than of order 1%. In that case, we 
can simplify equations (60)-(62) above by setting fa = and 
7^ = 1, to yield 

= n-Ari;, (64) 



where 



c dt fanH, P c + F 7 

N 1 (tR)-^r i I a- i n H 2Ca B Xefi/?> 
47Tfl 2 r? 



(65) 



In terms of the comoving volume , V/ = Airr] /3, equations (64) 
and (65) yield 

dVi _ A r 7 (ffi)-Ca B XeffnH, i a" 3 V/ 

dt f3in H ,i+ j^t- c [N 7 (t R )-Ca B X e m 2 Hti a- 3 V,] ' 



(66) 
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Equation (66) can also be derived from the cosmological equiv- 
alent of the photon conservation equation (23). This last ap- 
proach was first proposed (in a simplified form that did not ac- 
count for recombinations or the presence of helium) by White 
et al. (2003) (in their Appendix). Equation (66) was previously 
derived, including the effect of recombinations, but not of he- 
lium, by Yu & Lu (2005) using a different approach from ours. 

Following the approach in Shapiro & Giroux (1987), we de- 
fine the dimensionless volume as 

}' = Vi/Vsj, (67) 

where Vsj is the initial Stromgren volume, Vs.i = 4nrg ( /3 = 

N~,,i( a BCxeff n H,i)~ l , an d the dimensionless time as 

x = t/ti, (68) 

which is scaled in terms of the initial time rather than in terms of 
the recombination time as w was above. Substituting these rela- 
tions into equation (66), and using the fact that, for the matter- 
dominated era, x = a 3 / 2 , we obtain (for $ = Xeff) 

dy = X[l(x R )-y/x 2 ] 

dx 1 + (qx-hP)(x/y) 2/ V(x R ) -y/x 2 ] ' 

where 

A = — — = tiCa B n Hti , (70) 

Xeff'reca 

is the ratio of the age of the universe to the recombination time 
when the source turned on, divided by XeS, 

rs,i r s .iXeffX 

q = - — = — : — , (71) 

£'rec,i Cli 

is the light-crossing time of the initial Stromgren radius, in units 
of the initial recombination time, and 

l(x) = N^/N^n) = N y (x)/N^ (72) 

is the ionizing photon luminosity in units of its initial value (i.e. 
l(x) = 1 for a steady source). For q = (c — > oo), equation (69) 
reduces to the nonrelativistic form derived by Shapiro & Giroux 
(1987), for which they gave an exact analytical solution for the 
case of a steady source, as follows: 

y m (x) = Xe x/x [xEi(2, A/x)-Ei(2, A)] , (73) 

where Ei(2,x) = ^p-dt is the Exponential integral of second 
order. This exact solution also applies to a flat universe with a 
cosmological constant in the matter-dominated era (Iliev et al. 
2005a). As described in § 3.1, we can use this exact solution 
of the NR I-front equations to derive the exact solution of the 
relativistic I-front equations (i.e. arbitrary q / 0) directly, if we 
replace the time t in the NR solution by the retarded time t R . 
In that case, the exact, analytical, relativistic I-front solution is 
given by the following implicit equations: 

y (x) = y NR (x R ) = Ae A /-« [ XR Ei(2, X/x R ) - Ei(2, A)] , (74) 

and 

x = x R +x 2 ' 3 y\ll(x R ) (-^-] , (75) 

\ A Xeff/ 

where we have used the fact that a(x) = x 2 / 3 , x = x R + n tP (x) / (ctf), 
and r/ iP (x) = a(x)r^,[ j yNR(xR)] 1,/3 . Equations (74) and (75) easily 
reduce to a cubic equation for x(x R ), which has an exact analyt- 
ical solution with one physical root. Since it is a rather lengthy 
expression, we omit it here. 



For each value of A, we obtain a family of solutions for dif- 
ferent values of q. In fiducial units 

, = 0.0217^^(1+^ (fig) " 5 , (77) 

where (1 +z,-)io = (1 +z,)/10, /V 7 , 56 = ^/(lO^s -1 ), and Sl b is 
the mean baryonic density in units of the critical density. The 

-2/3 

solutions for A = 0.5, 1,4, and 10 and qx e g = 0,0.1,0.2 and 
0.5 are shown in Figures 5 and 6. The relativistic correction 
is most significant for small A and large q values, i.e. at later 
times, stronger sources and little or no clumping of the gas. 
This leads to somewhat counterintuitive behavior. For example, 
increasing the clumping factor increases q, and naively, based 
on equation (69), could be expected to cause a stronger rela- 
tivistic correction. However, it also increases the value of A, 
which results in a smaller relativistic correction. This is be- 
cause larger clumping of the gas increases the recombinations 
in the H II region, which decrease the ionizing flux at the front, 
slowing it down into the nonrelativistic regime. Below we shall 
see that infall around the source has a similar effect, increas- 
ing the density near the source and rendering the relativistic 
corrections somewhat weaker than a simple, uniform-density 
treatment would predict. 

As in the case of a spherical I-front in a static medium dis- 
cussed in §3, the apparent size of a cosmological H II region 
observed at a given time will depend upon the angle between 
the line of sight to the source and the line connecting the source 
to the point we observe on the I-front surface. This difference 
results from the difference in the light travel times for light re- 
ceived by the observer from different locations on the front. If 
we were to measure the radius of the H II region transverse 
to the line of sight at the distance of the source, for example, 
and interpret this as the proper radius given by the relativis- 
tic solution at some time t, then the line-of-sight radius on the 
near-side of the source would be given by the nonrelativistic 
solution at the same physical time. Hence, the line-of-sight ra- 
dius will appear to be larger than the transverse radius and will 
appear to have expanded toward the observer at superluminal 
speeds. This deviation of the shape of the cosmological H II 
region from spherical was also discussed by Yu (2005) in the 
context of the H II regions around QSO's at redshifts z > 6 
whose spectra show Gunn-Peterson troughs. 

5. COSMOLOGICAL RELATIVISTIC I-FRONTS IN PRESENCE 
OF SMALL-SCALE STRUCTURE AND INFALL 

When self-shielded neutral structures (minihalos) are present 
on the path of the propagating cosmological I-front, we fol- 
low the approach in Iliev et al. (2005a). In this case the I-front 
jump condition in equation (60) and the I-front evolution equa- 
tion (66) are modified by setting 

& = ( 1 - /coiOXeff + [ 1 +A(He)]/ coll ,MH£ (78) 
where / co u is the total collapsed baryon fraction, / co ii,mh is the 
collapsed fraction of just the minihalos, and £ is the number 
of ionizing photons consumed per minihalo atom such that £ is 
the appropriate average over the distribution of minihalos at the 
instantaneous location of the global I-front. This photon con- 
sumption factor is flux-dependent and in the relativistic case the 
appropriate flux is F, the flux of ionizing photons in the frame 
of the source, as the relative velocity between the source and 
the minihalos is much smaller than c. 
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FIG. 5. — Cosmological I-front radius versus time for steady source in the mean IGM, with comoving radius in units of the time- varying, comoving Stromgren 
radius, rs(t) = a{t)rs(U), and time in units of the age of universe f, at source turn-on. Different curves in each panel correspond to different values of q = rsj/(ct re c,i), 
with g = (i.e. nonrelativistic limit, top line) and (from top to bottom) q = 0.1,0.2 and 0.5. Each panel is for different values of A = f;/(Xeft*rec.i), as labeled. We 
have assumed that Xeff = 1 for simplicity. 




FIG. 6. — Same as Figure 5 but for the I-front peculiar velocity, instead, where vi, pec is in units of (r si /ti). The speed of light, c, corresponds in these units 
(XeffX/q). We assume Xeff = 1 f° r simplicity. 
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FIG. 7. — Cosmological I-front radius in present (z = 0) comoving units, r c = r/(l +z,), versus time (since source turn-on) for a large starburst (left panels) and a 
luminous QSO (right panels) in the local IGM with infall density profile, both starting at Zi = 7, with assumed parameters as described in the text, for the cases with 
mean (C = 1) and clumped (C = 10) IGM, with and without minihalos present, as labeled, for the relativistic case (solid) and the corresponding nonrelativistic one 
(dashed). 



The average number of ionizing photons absorbed per mini- 
halo atom at the current location of the I-front, £, is given 
by integrating the photon consumption rate £ per individual 
minihalo over the Press-Schechter mass function of miniha- 
los at that epoch. The photon consumption rates for individual 
minihalos were determined by Iliev et al. (2005b) using high- 
resolution, numerical gas dynamics simulations that included 
radiative transfer. They have summarized their detailed results 
by providing fitting formulae for the individual minihalo con- 
sumption rates and their dependence on halo mass, redshift, in- 
cident flux and ionizing source spectrum. According to Iliev 
et al. (2005b), 

&M,Z,F ) = 1 + 4. 4M ^34 + 0.0231o glo M 7jF 0.199-0.0421o gl oF 0(1 ^ 

(79) 

for a 5 x IQ 4 K blackbody external ionizing spectrum, represent- 
ing Pop. II stars (and similar numbers for a QSO spectrum), for 
a minihalo of mass M-j (in units of 10 7 M Q ), overtaken by an 
intergalactic I-front at a redshift of z which is driven by an ex- 
ternal source of flux Fq, the flux in units of that from a source 
emitting 10 56 ionizing photons per second at a proper distance 
d of 1 Mpc, i.e. Fq = N^^/dl^ . Again, the relevant flux to 
use in calculating £ is F, the photon flux in the source frame, as 
it is only the ionization front, not the self-shielded structures, 
that is moving relativistically with respect to the source. Our 
model takes account of the biasing of minihalos with respect to 
the ionizing sources (i.e. the clustering of minihalos in space) 
by modifying the mass function of the minihalos as a function 
of their distance from the central source halo in each H II re- 
gion. It also includes a simple treatment of the effect of local 
infall on the density of the IGM and the clustering of minihalos 
surrounding each source [see Iliev et al. (2005a) for details]. 

6 Our I-front propagation equation remains valid even if the source luminosity decays, as long as the I-front speed remains positive and the luminosity decay time 
is greater than the ionization equilibration time inside the H II region. This equilibration time is much smaller than the fully-ionized recombination time, however, 
since it is given by the product of that recombination time and the small neutral fraction inside the H II region. The propagation equation is also valid even when it 
yields a negative I-front velocity, but only if the luminosity decay time exceeds the fully-ionized recombination time. In that case, the I-front shrinks to match the 
instantaneous Stromgren sphere. 



In the results below we specialize to a Cold Dark Matter 
(CDM) cosmological model with parameters h = 0.7, flo = 0.3, 
J7a = 0.7, and = 0.045, erg = 0.87, where 51a is the total vac- 
uum density in units of the critical density, a\ is the variance 
of linear fluctuations extrapolated to the present and filtered on 
the 8/z _1 Mpc scale, and n p is the index of the power spectrum 
of primordial density fluctuations (Spergel et al. 2003). Note, 
however, that as the fronts we are considering are all at high red- 
shifts, the value of fl\ has no direct impact on our calculations. 
Finally, the Eisenstein & Hu (1999) transfer function was used. 
In Figure 7, we show our results for the I-front radius evolu- 
tion for two representative sources. The first is a galaxy of total 
mass lO n M at an initial redshift of z, = 7 (corresponding to a 
3-er density fluctuation), which hosts a starburst that produces 
a total of / 7 = 250 ionizing photons per baryon within a life- 
time of t s = 3 Myr, corresponding to N 7 = 4.08 x 10 55 s _1 . For 
t > t s , the ionizing photon production is assumed to drop as f " 4 5 
(Haiman & Holder 2003). The second source roughly corre- 
sponds to a high-redshift luminous QSO, similar to those found 
by the Sloan Digital Sky Survey (e.g. Fan et al. 2004). Here 
we assume an ionizing photon production rate of N 7 = 10 57 s _1 , 
a source lifetime t s = 30 Myr starting at z = 7, and a host 
halo of 10 13 M Q . For t > t s , we model the drop in the ioniz- 
ing photon production rate as an even steeper power of time, 
N 7 (t > t s ) oc f~ 10 , although the precise form of this fall-off is 
unimportant since we are interested in the expansion phase of 
the H II region. 6 We show the relativistic results (solid) and 
the naive nonrelativistic results (dashed) for direct comparison, 
and consider four cases with mean (C = 1) and clumped (C = 10) 
IGM, with and without minihalos present. 
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For both sources, the I-front expansion is initially superlumi- 
nal in the nonrelativistic model, causing the H II region radius 
to be significantly larger than the correct, relativistic model. 
These differences diminish somewhat at later times, but remain 
significant throughout the expansion phase of the H II region. 
For both sources, the IGM clumping has only a modest effect, 
since even for C = 10 the recombination times are longer than 
the source lifetimes. The presence of minihalos, on the other 
hand, decreases the size of the H II regions more considerably, 
by about 25% in each case. We have also calculated the I-front 
expansion for a weaker source at higher redshift (a galaxy halo 
of total mass 10 8 M Q at initial redshift z, = 15, corresponding to 
a 3 — a density fluctuation, which hosts a Pop. II starburst that 
produces a total of / 7 = 250 ionizing photons per baryon within 
a lifetime of t s = 3 Myr) and, as expected from Figure 1, the 
relativistic correction in this case was negligible. 

Finally, we quantified the impact of relativistic effects on the 
global progress of reionization by recalculating the reionization 
models described in detail in Iliev et al. (2005a). For all mod- 
els, the relativistic corrections were small, decreasing the global 
ionized mass fraction at any given time by ~ 1 % as compared 
to the nonrelativistic treatment. The corresponding effect on 
the mean electron-scattering optical depth r was even smaller, 
never exceeding a small fraction of 1%. Thus while the rela- 
tivistic expansion of an I-front about an individual strong source 
close to overlap is important, it can generally be ignored when 
computing the global reionization history, especially given the 
other significant uncertainties involved. 

6. SUMMARY AND CONCLUSIONS 

We have considered the effect of special relativity on the 
propagation of ionization fronts. The standard (nonrelativis- 
tic) treatment of I-fronts neglects the finite speed of light and 
the relative motion of the radiation source and the gas it ion- 
izes. As a result, when the ionizing photon-to-atom ratio ex- 
ceeds unity at the location of the I-front, unphysical superlumi- 
nal I-front motion is predicted to occur. We have generalized 
the I-front continuity jump condition and the equations that de- 
scribe the transfer of ionizing radiation between the source and 
the I-front to take account of these things in a relativistically- 
correct way. By combining these equations, we have derived 
the relativistic equation of motion for the R-type I-front created 
when a point source of ionizing radiation turns on in a neutral 
gas, including the possibility that the gas is non-uniform and 
in motion. We then applied this relativistic equation to several 
illustrative problems of interest in the theory of H II regions 
in the interstellar and intergalactic media, including the cos- 
mological problem of the reionization of the universe. These 
include the cases of a steady source in a static gas that is ei- 
ther uniform, spherically-symmetric with a power-law density 
profile or plane-stratified. We have also solved the problem of 
a point source in a cosmologically-expanding IGM. We have 
shown that, in all cases, the radius of the relativistic I-front at 
a given physical time is identical to that in the nonrelativistic 
solution in which the physical time t is replaced by the retarded 
time, equal to t minus the light-travel time from the source to 
the front. 

In the case of a static medium, we obtain an exact analytical 
solution for the problem of a steady source in a uniform gas, 
which depends only on the parameter q, the time for light to 
cross the Stromgren radius, in units of the recombination time. 
Applying this solution to the problem of an O-type star or OB 



association in a molecular cloud, we find q values of order unity 
over the full range of appropriate densities and fluxes. How- 
ever, while this implies that corrections are formally important, 
they only affect ionization fronts during the rapidly-expanding 
R-type phase, which lasts only until the I-front expands to fill 
the Stromgren sphere, and not during the slower D-type hydro- 
dynamic expansion that dominates the bulk of the evolution of 
these H II regions. 

For a source in a nonuniform, static gas, however, I-fronts 
can, in principle, accelerate back up to relativistic speeds if 
they propagate down a steep enough density gradient. We have 
solved this problem for the cases of a power-law density profile 
in spherical symmetry and a plane-stratified gas. These solu- 
tions depend both upon the value of the parameter q calculated 
at the central density in the core of these nonuniform profiles 
and on the dimensionless ratio of the core radius of the density 
profile (or the scale-height in the plane-stratified case) to the 
Stromgren radius. Such solutions are relevant, for example, to 
the "champagne phase" of the H II region, which results when 
the D-type I-front from a stellar source in a dense environment 
detaches from the shock that leads it during the dynamical ex- 
pansion phase, and the resulting R-type front runs away towards 
large radii. In that case, relativistic corrections are important, 
not only right after the source turns on, but at late times as well. 

In the cosmological case, we have derived the exact, analyt- 
ical solution for the dependence of the relativistic I-front ra- 
dius on time, by generalizing the solution derived by Shapiro & 
Giroux (1987). We find a complete family of solutions that are 
described by two dimensionless parameters: q, which is now 
defined as the light crossing time at the Stromgren radius in 
units of the initial recombination time, and A, the age of the 
universe at source turn-on in units of the initial recombination 
time. In this case, relativistic effects are again important only if 
q :> 0.1 , but now only if A <; 1 as well. Thus relativistic correc- 
tions are significant only in cases in which both the source flux 
is large (corresponding to a large q) and the redshift is relatively 
low (corresponding to a small A). 

Incorporating relativistic effects into the more exact model- 
ing described in Iliev et al. (2005a) allows us to assess the im- 
portance of relativistic effects in a realistic reionization scenario 
that includes both infall and the impact of small-scale (self- 
shielding) structures. In the individual case of a large galaxy or 
luminous QSO at the end of the reionization epoch (z ~ 7), rel- 
ativistic effects can be quite significant: resulting in m 4 and 2 
times smaller ionized volumes after 1 and 3 Myrs, respectively, 
in the galactic case, and w 20 and 2 times smaller volumes af- 
ter 1 and 10 Myrs, respectively, in the QSO case. The relative 
differences then decrease with time, but remain significant until 
the source luminosity decays. Afterwards, at time t ~ 2t s , the 
relativistic solution for the H II region radius finally matches 
the nonrelativistic one. Although the I-front propagation be- 
comes subluminal earlier than this, and thus the usual nonrela- 
tivistic expression for its speed becomes correct, the radius of 
the H II region remains smaller throughout the source lifetime. 
Furthermore, the presence of small-scale structure can decrease 
the size of the H II region by up to an additional ^25%. 

Generalizing these results to the issue of global reionization, 
however, yields far less dramatic changes. As most of the inter- 
galactic medium was ionized by smaller (lower q) and higher- 
redshift (higher A) sources, we find that relativistic corrections 
change the ionized mass fraction at any give time by ~ 1 % with 
respect to a nonrelativistic approach, and relativistic corrections 
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to the mean electron-scattering optical depth are even smaller. 
Thus I-fronts, while an intrinsically relativistic phenomenon, 
can be accurately treated nonrelativistically in studies of global 
cosmological reionization, but not for the isolated H II regions 
of short-lived luminous sources at late times, close to the epoch 
of overlap. 

Finally, we note that the original nonrelativistic, cosmolog- 
ical I-front solutions of Shapiro & Giroux (1987) correctly re- 
late the observed line-of-sight distance of the I-front from the 
source which created it to the source luminosity and IGM neu- 
tral density, even for cases in which those nonrelativistic so- 
lutions yield a superluminal expansion speed. Thus, the inter- 
pretation of the line-of-sight apparent sizes of the cosmological 
H II regions recently observed in the IGM surrounding high- 
redshift sources like the SDSS quasars at z > 6 does not re- 
quire the relativistic correction described here. This somewhat 
surprising result was noted previously by White et al. (2003) 
and others (e.g. Yu & Lu 2005). Our results help to illuminate 
the connection between the nonrelativistic and relativistic treat- 



ments, however, and distinguish the circumstances under which 
the relativistic corrections do make a difference. For example, 
when we observe a distant source and match the apparent size 
of its H II region along the line of sight to the nonrelativistic 
solution evaluated at some time t (as measured from the source 
turn-on time), the apparent transverse size at the distance of the 
source will then be given by the relativistic solution evaluated 
at this same time t . 



ITI would like to express his sincere thanks for the hospital- 
ity shown to him by the Kavli Institute for Theoretical Physics, 
during part of this research. ES was supported by the National 
Science Foundation under grant PHY99-07949. This work 
was partially supported by NASA Astrophysical Theory Pro- 
gram grants NAG5- 10825 and NNG04G177G and Texas Ad- 
vanced Research Program grant 3658-0624-1999 to PRS. MAA 
is grateful for the support of a Department of Energy Graduate 
Fellowship in Computational Sciences. 



REFERENCES 



Bruscoli, M, Ferrara, A., & Scannapieco, E. 2002, MNRAS, 330, L43 
Churchwell, E. 2002, ARA&A, 40, 27 

Ciardi, B., Ferrara, A., Governato, F, & Jenkins, A. 2000, MNRAS, 314, 61 1 
Ciardi, B., Stoehr, F, & White, S. D. M. 2003, MNRAS, 343, 1101 
Eisenstein, D. J. & Hu, W. 1999, ApJ, 511, 5 
Evans, N. J. 1999, ARA&A, 37, 311 
Fan, X., et al. 2004, AJ, 128, 515 

Franco, J., Tenorio-Tagle, G., & Bodenheimer, R 1990, ApJ, 349, 126 

Furlanetto, S. R., Zaldarriaga, M., & Hernquist, L. 2004, ApJ, 613, 1 

Garay, G. & Lizano, S. 1999, PASP, 111, 1049 

Gnedin, N. Y. 2000, ApJ, 535, 530 

Haiman, Z. 2002, ApJ, 576, LI 

Haiman, Z. & Holder, G. P. 2003, ApJ, 595, 1 

Hosokawa, T. & Inutsuka, S. 2005, ApJ, 623, 917 

Iliev, I. T., Scannapieco, E., & Shapiro, P. R. 2005a, ApJ, 624, 491 

Iliev, I. T. & Shapiro, P. R. 2001, MNRAS, 325, 468 

Iliev, I. T., Shapiro, P. R., & Raga, A. C. 2005b, MNRAS, 361, 405 

Kogut, A., et al. 2003, ApJS, 148, 161 

Lebron, M. E., Rodriguez, L. F, & Lizano, S. 2001, ApJ, 560, 806 



Lee, Ft., Chen, W. P., Zhang, Z., & Hu, J. 2005, ApJ, 624, 808 

Lefloch, B., Lazareff, B., & Castets, A. 1997, A&A, 324, 249 

Malhotra, S. & Rhoads, J. E. 2004, ApJ, 617, L5 

Mathews, W. G. & O'Dell, C. R. 1969, ARA&A, 7, 67 

Oh, S. P. & Furlanetto, S. R. 2005, ApJ, 620, L9 

Ricotti, M., Gnedin, N. Y, & Shull, J. M. 2002, ApJ, 575, 49 

Santos, M. R. 2004, MNRAS, 349, 1137 

Shapiro, P. R. & Giroux, M. L. 1987, ApJ, 321, L107 

Shapiro, P. R., Iliev, I. T., & Raga, A. C. 1999, MNRAS, 307, 203 

— . 2004, MNRAS, 348, 753 

Shapiro, P. R., Milgrom, M., & Rees, M. J. 1986, ApJS, 60, 393 

Sokasian, A., Abel, T., Hernquist, L., & Springel, V. 2003, MNRAS, 344, 607 

Spergel, D. N, et al. 2003, ApJS, 148, 175 

Stern, D., et al. 2005, ApJ, 619, 12 

White, R. L., Becker, R. H., Fan, X., & Strauss, M. A. 2003, AJ, 126, 1 

Wyithe, J. S. B. & Loeb, A. 2003, ApJ, 586, 693 

— . 2004, Nature, 432, 194 

Yu, Q. 2005, ApJ, 623, 683 

Yu, Q. & Lu, Y. 2005, ApJ, 620, 31 



